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Abstract: Multilayered microresonators commonly use sensitive coating or piezoelectric 
layers for detection of mass and gas. Most of these microresonators have a variable 
cross-section that complicates the prediction of their fundamental resonant frequency 
(generally of the bending mode) through conventional analytical models. In this paper, we 
present an analytical model to estimate the first resonant frequency and deflection curve of 
single-clamped multilayered microresonators with variable cross-section. The analytical 
model is obtained using the Rayleigh and Macaulay methods, as well as the Euler-Bernoulli 
beam theory. Our model is applied to two multilayered microresonators with piezoelectric 
excitation reported in the literature. Both microresonators are composed by layers of seven 
different materials. The results of our analytical model agree very well with those obtained 
from finite element models (FEMs) and experimental data. Our analytical model can be 



Sensors 2011, 11 



8204 



used to determine the suitable dimensions of the microresonator's layers in order to obtain a 
microresonator that operates at a resonant frequency necessary for a particular application. 

Keywords: bending resonant frequency; Euler-Bernoulli beam theory; Macaulay method; 
multilayered microresonator; Rayleigh's method 



1. Introduction 

Multilayered microresonators fabricated using microelectromechanical systems (MEMS) have 
potential applications such as mixer- filters [1,2] and detection of mass [3,4], prostate-specific 
antigens [5], proteins [6], virus [7,8], chemical species [9], gas [10,11], and magnetic fields [12-14]. 
These microresonators have several advantages such as small size, high sensitivity, high resolution, 
low power consumption, and low cost in mass production [15,16]. They can be excited in the first 
bending resonant frequency through different types of sources such as magnetic, electrostatic, thermal, 
and piezoelectric [17]. Their frequency shift and deflections can also be monitored using capacitive, 
piezoelectric, piezoresistive or optical sensing techniques [18,19]. 

The microresonators are extremely sensitive to surface processes due to their large surface area to 
mass ratios [20]. For chemical and biological sensing, the surface of microresonators can be covered 
with a sensitive coating to absorb the chemical and biological species [21,22]. Several microresonators 
also use piezoelectric multilayers due to their advantages such as compact structure, low driving 
voltage, self-actuation and self-sensing capability, and easy integration with electronic circuits [23]. 

An important issue in the mechanical design process of a multilayered microresonator with variable 
cross-section is the correct determination of its fundamental resonant frequency. Another important 
parameter is the estimation of its deflection under an excitation load, which can be used to study its 
mechanical behavior. Theoretical models can predict the resonant frequency shift of microresonators 
caused by variations of their geometrical variables and mechanical properties. Lobontiu etal. 
developed lumped and distributed-parameter models for predicting the bending or torsion resonant 
frequencies of microresonators with variable cross section and serially connected [24-28]. They also 
obtained analytically the deflection of this microresonator type; nevertheless, their models do not 
consider multilayered microresonators. Dufour and Fadel [29] reported analytical expressions in order 
to calculate the resonant frequencies of thin parallelepiped microcantilevers formed by a structural 
layer and a sensitive coating. They calculated the bending resonant frequency of microcantilevers 
using the deflection expression of its bending resonance mode. Their models are suited to 
microcantilevers with a rectangular plate and a sensitive coating. Ferguson et al. [30] presented an 
analytical model for predicting the fundamental bending resonant frequency of a composite free-free 
microbeam, which is connected to four torsional anchors. In addition, they determined the fundamental 
torsional frequency of the composite anchors. However, these models are limited to composite 
free-free microbeams and anchors with uniform cross section. Delia et al. [31] employed an analytical 
model in order to study the free vibration of de laminated multilayered beams with two fixed ends. 
They analyzed the fundamental frequency shift of clamped-clamped multilayer beams with single and 
double delaminations in terms of the following parameters: normalized axial stiffness, normalized 
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bending stiffness, and the relative slenderness ratio. This model is limited to delaminated multilayered 
beams with uniform width. Sampath et al. [32] determined the resonant frequency of a two-layer 
hybrid microcantilever with uniform cross section, which is integrated by an elastic base and 
viscoelastic layer. They reported that the molecular absorption process affects the resonant frequency 
of the microcantilever due to that modifies the mass load and the viscoelastic properties of the 
sensitive coating. Li et al. [33,34] studied the nonlinear transverse vibrations of a clamped-clamped 
composite microbeam using a refined integro-partial differential model. This model considers a 
composite microbeam with stepwise axially varying properties, rectangular cross section, and non-flat 
equilibrium position {i.e., a buckled beam). Based on this model, they estimated the first natural 
frequency of a clamped-clamped composite microbeam. For this, they employed a buckling factor that 
depends of the static compressive forces, stiffness, and geometric properties of the microbeam. Finally, 
they only reported the natural frequency of a multilayer microbeam with two clamped ends, which has 
a uniform width and a non-flat equilibrium position. Later, Li et al. [35] adopted a curved cross-section 
model to refine their previous buckling microbeam model [33,34]. This new model considers the 
cross-section curvature of the clamped-clamped composite microbeam induced by residual stress, 
which affects its structural stiffness and spatial response. Using this model, the analytical resonant 
frequencies of the clamped-clamped composite microbeam is closer to the experimental observations. 
However, it cannot be used for multilayered microresonators with variable width or serially connected 
multisegments. Pasini [36] developed a bending response model for multilayered microresonators 
considering their cross-section shape, symmetry and layers number, as well as their materials 
properties. However, his model does not consider different microbeam types connected in series. 
Edqvist et al. [37] established a general theoretical model for studying the quasi-static and dynamic 
electromechanical responses of piezoelectric multilayered microbeams with one fixed end. Their 
model was obtained through Euler-Bernoulli theory and only considers layers with uniform cross 
section. Herrera-May et al. [38,39] used theoretical models to estimate the first bending resonant 
frequency of microresonators integrated by an arrangement of microbeams with different 
cross-sections. All these models have important characteristics that can be useful in the design phase of 
some microresonators. However, several they have limitations to predict the resonant frequency of 
multilayered microresonators with variable cross-section. In order to overcome all these drawbacks, 
we present an analytical model for the first bending resonant frequency of multilayered 
microresonators with variable cross section and one fixed end. This model includes microresonators 
composed by an arrangement of different multilayered microbeam types. Furthermore, our model can 
be used to estimate the deflection expression of the proposed microresonators. The analytical model is 
based on the Rayleigh's method, Macaulay's functions, and Euler-Bernoulli beam theory. It does not 
consider the damping, residual stress and shear deformation of the microresonators. 

This paper is organized as follows: in Section 2, we describe the analytical model for the bending 
resonant frequency of a multilayered microresonator with variable cross section. Next, in Section 3, 
this model is applied to two multilayered microresonators with mass sensing applications reported in 
the literature. The results of our analytical model are compared with those obtained from a simple 
cantilever model, FEMs, and experimental data. Finally, we presented our conclusions and proposed 
future research in Section 4. 
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2. Analytical Modeling 

This section describes the analytical modeling to predict the first bending resonant frequency and 
deflection of multilayered microresonators with variable cross section. Rayleigh's method can predict 
the first bending resonant frequency of structures with complex geometrical shapes considering energy 
conservation in the structures [40]. This method requires the deflection expression, bending stiffness, 
and mass per unit length of the structures. 

The bending deflection response y(x,t) at a given point of a resonant structure with a harmonic 
motion can be determined as a product between a spatial function y(x) and a time-dependent one [27]. 
Thus, y(x,t) = y(x)sm(2%ft), where t is the time and /is the frequency. 

For out-of plane bending in single and double clamped beams with length L and cross-section area 
A(x), the maximum potential energy (U m ) and kinetic energy (K m ) are determined by [41]: 



where E, I, and p are the Young's modulus, inertia moment, and density of the beam, respectively. 

Based on Rayleigh's method, the first bending resonant frequency (f r ) of single and double clamped 
beams can be obtained assuming energy conservation (U m = K m ) in the beams. By substituting 
Equations (1) and (2) into U m = K m , the f r is given by: 



In order to use the Rayleigh's method in a multilayered microresonator with variable cross section, 
we need to know its elastic centroid, bending stiffness and mass per unit length. We propose a 
multilayered microresonator with variable cross section, as shown in Figure 1 . This microresonator is 
formed by layers of different materials and geometric shapes, which are symmetrically distributed with 
respect to the xy-plane. The microresonator is divided into four different sections, which contain mth, 
nth, pth., and qth layers, respectively. We employed the nomenclature y'th section in order to describe 
each one of the four sections of the microresonator. Figure 2 shows the geometrical nomenclature 
proposed for Mi layer placed over any of the y'th section. In this case, b t sj, and % represent the width 
and thickness of the z'th layer located on the y'th section. In addition, the parameter hisj indicates the 
distance from the bottom plane of the first layer to the top plane of the z'th layer placed over the y'th 
section. Based on our multilayered microresonator, several microresonator designs used for detection 
of mass and gas can be obtained. Figure 3(a-d) illustrates the configurations of four different 
multilayered microresonators that are derived from our general-microresonator design. An analytical 
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model for estimating the first bending resonant frequency of our general microresonator can help 
designers and researchers in the design phase of particular microresonators. For instance, an analytical 
model for a general multilayered microresonator can be used to know the effect of the materials and 
dimensions of the multilayers on the resonant frequency and deflection of particular microresonators. 

Figure 1. View of the multilayered microresonator with variable cross section proposed in 
this work. 




Figure 2. Geometrical nomenclature proposed for the kth layer located on the y'th section of 
the proposed multilayered microresonator. 
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Figure 3. Examples of multilayered microresonators obtained from the proposed 
multilayered microresonator. These microresonators could are composed by silicon layers 
with (a-b) sensitive coating and (c-d) piezoelectric layers. 




Figure 4 shows the reaction loads, bending moments, and uniformly distributed loads of the 
proposed microresonator. 

Figure 4. View of the load types acting on the proposed multilayered microresonator. 




For this case, Ro and Mo represent the total reaction load and bending moment at the fixed end of 
the microresonator. Furthermore, CDsj indicates the uniformly distributed weight of all the layers located 
in the y'th section. The analytical model developed in this work is based on the following assumptions: 

i) Multilayered microresonator is integrated by layers made of homogeneous and isotropic 
materials. The film mechanical properties of the different layers must be considered; 

ii) The microresonator layers are symmetrically distributed with respect to xy-plane; 



Sensors 2011, 11 



8209 



iii) The bending vibration of the microresonator occurs along of the j-axis and its vibration 
amplitude is smaller than its total thickness; 

iv) The plane sections of layers do not deform, i.e., transverse shear strain is neglected; 

v) The xz-plane of the coordinate system is located at the elastic centroide of the microresonator; 

vi) The total length of the microresonator is greater than 10 times its total thickness. This is a 
necessary condition to use the Euler-Bernoulli beam theory; 

vii) Damping and surface effects (e.g., surface energy, surface tension, surface relaxation, and 
surface reconstruction) are neglected; 

viii) The residual stress in the microresonator is neglected; 

ix) The nonlineality and viscosity of the layers are not considered; 

x) The undercut to the microresonator caused by etching processes is neglected. 

For our microresonator, the elastic centroid (as,) in the y'th section of the microresonator can be 
determined as [42]: 



where t iSj = htsj ~ hy-iyy, the subscript 5} is related to the y'th section, the value k depends of the layers 
number located over each one of the four sections (i.e., k can take either values m, n, p 
or q), A s symbolizes the domain occupied by the y'th section, ysj(x) is the deflection along of the y'th 
section, Etsj is the Young's modulus of the z'th layer located in the y'th section, /z(,-i)# is distance from the 
bottom plane of the first layer to the top plane of the (z'-l)th layer placed over the y'th section, (ES) S is 
the elastic modulus weighted first moment of area and {EA) S the elastic modulus weighted area. For 
the first section, bisi considers the overall width of the z'th layer (i.e., btsi includes the sum of the width 
of the three layers located on the same distance h t si, which have the same material). In addition, the 
term hosj = 0. 

The effective bending stiffness (EI z )sj in the y'th section of the multilayered microresonator can be 
calculated as: 




(4) 
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The maximum potential energy (U m L) and kinetic energy (K mL ) of the multilayered microresonator 
can be obtained as: 




K 



mL 



if m If 2 if " |f 2 

v !=i / o v !=i / i Si 

if p XT if ? VT 



(6) 



(V) 



where L S u = L S \ + L S i, L S m = L S i + L S2 + L S3 , L si2 i4 = L S \ + L S2 + L S3 + L s4 , and X = 2nf. 

Using the Rayleigh's method, the first bending resonant frequency of the multilayered 
microresonator is obtained as: 



fr = 



1 

In 



U, 



,1/2 



mL 



\ K mL< ' ^ ) 



(8) 



Equations (6-8) show a dependence of the bending resonant frequency with respect to the 
deflection y(x), mechanical properties (Young's modulus and density) and geometrical parameters of 
the layers from each section of the resonant microestructure. Any variation of these factors will cause a 
change in its resonant frequency. For the case of mass and gas sensing, the particles deposited on the 
surface of a microresonator can affect its bending resonant frequency due to changes of its mass and 
stiffness. This frequency shift is used to detect the amount and type of particles deposited on the 
microresonator surface. For the case of small particles deposited on a concentrated region of the 
microresonator surface, the resonant frequency shift can be estimated only considering the kinetic 
energy variation. In order to determine this energy shift, the mass added by the particles must be 
included into Equation (7). For instance, if the mass of the particles is considered as a small 
concentrated mass load im a dd) located to a distance x a dd of the microresonator's fixed end, then the 
term m a ddy (x a dd)/2 must be summed into Equation (7). For more information of this topic we 
recommended the work developed by Dong et al. [11]. We do not include the study of the resonant 
frequency shift of the microresonator. 

The deflection ysj(x) in each one of the four sections of the microresonator are obtained using 
Euler-Bernoulli beam theory and Macaulay's method [43]. This method is suitable to describe several 
load types acting on structures with variable cross section [44]. For this, it employs functions whose 
nomenclature considers a bracket notation (x - c) £ , since these functions have the zero value for x < c 
and (x - c) e for x > c [45]. The value of superscript e depends of the load type acting on structure. For 
instance, e = 0 for a uniformly distributed load, £ = -1 for a concentrated load, and £ = —2 for a 
bending moment. 
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Thus, the deflection in each one of the four sections of proposed microresonator are determined 
from: 

, d 2 Vo (x) 

( E h) Si ^^ l = M Si {x) 0<x<L S] (9a) 
(EI z ) s d -^-=M Si {x) L Si <x<L Sii (9b) 

where M$j is the bending moment in the jth section of the multilayered microresonator, which can be 
determined through Macaulay's functions. 

The deflection ys/x) of our microresonator must satisfy the following boundary conditions: 

A(0) = 0 

y Sl ( L s 1 )=ys 2 ( L s l ) 
y S2 ( L sJ=ysS L s n ) 

3^5 3 C^x 123 ) ys 4 ^^S m ) 

Using Macaulay's method, we find the total load function F(x) of the microsensor. It can be 
described as [43]: 

F(x)=-M 0 (x-0}~ 2 +R 0 {x-0}~ 1 -0) Sj (x-0)° +co Sj {x-L^ )° -co Si (x-L^ )° +co Si (x-L^ }° 

/ \° / \° / \° / \° 

-co, (x-L, ) +co. (x-L~ ) -co, (x-L, ) + co, (x-L, ) 

°3 \ J 12 / J 3 \ J 123 / J 4 \ ^122 I J 4 \ J 1234 / 

with 

4 

^=2>*A ( 12 ) 



3^(0) 


= 0 


dx 




dy Sl ( L s) 




dx 


dx 






dx 


dx 


^ys 3 (^s 123 ) 




dx 


3x 



(10) 



1=1 



4 f 1 

0 S,- S, ^ 2 ^ S j-i V2 V: 

k 



(13) 



HP iSl Sb iS t iSi (14) 



7 . , j y j 
1=1 



where represents the weight per unit length in the y'th section and g is the gravitational acceleration. 
For the first section, bis\ considers the overall width of the all layers located in h t s\ (e.g., if there are 
three layers made with the same material located to a distance htsi, then bjsi includes the sum of each 
width of the layers). 
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The shear load function V(x) is obtained by integrating Equation (11) with respect to x and using the 
integration rules of the Macaulay's functions [43]. The function V{x) is given by: 

V(x) = -M 0 (x-O)' 1 +R, (x-0)° -0) Si (x-0) 1 +0) Si (x-L S] )' -co Si (x-L^ )' +co Si (x-L^ )' 

-fl\ (x-L<; ) +a)„ lx-L, ) -co, ix-L, ) + co<, Ix-L^ ) +C 

J 3 \ J 12 / J 3 \ J 123 / J 4 \ J 123 / J 4 \ J 1234 / 1 



Integrating Equation (15) with respect to x, the bending moment function M(x) is determined as: 

1 1 2 1 2 1 2 

M(x) = -M 0 (x-0}° +R o (x-0)' --co Si {x-0} 2 +-co Si (x-L Si ) --co Si (x-L s ) +-co Si (x-L Sn ) 

1 21 21 21 2 0^) 
— <y c (x-L, ) H — co, (x-L, ) — <u, (x-L<, ) h — co, (x-L, ) +Cx+C 

2 j 3 \ °i2 / 2 2 3 2 34 1 

where C\ and C2 are constants calculated through boundary conditions (V(0) = Ro and M(0) = Mo) of 
the shear load and bending moment functions. Substituting these conditions into Equations (15) 
and (16), we obtain C\ = C2 = 0. 

The bending moment functions for each one of the four sections of the multilayered microresonator 
are determined from Equation (16): 



For 0<x<L s 



M Si (x) = -M 0 (x-0)° +R 0 (x-0} 1 ~co Si {x-0f (i 7a ) 



For L s <x<L s 



12 



1 1 2 1 2 

M S2 (x)=-M 0 (x-0)° +R 0 (x-0) 1 --co Si {x-0) 2 +-co Si (x-L Si ) --co Si (x-L s ) (17b) 

For L, <x<L, 

°12 °123 

M S} (x) = -M 0 {x - 0)° + R, (x - 0)' - ]- co Si {x - of + 1 Q) S] (x - L s< | 



--co, (x-L,\ +— co, (x-L, ) --co, (x-L, v 

2 J2 \ / 2 2 \ 12 / 2 1 



(17c) 



For L„ <x<L? 

J 123 J 12 



1 1 2 1 

M 54 (x)=-M 0 (x-0) 0 +^(x-0) 1 --^ | (x-0> 2 +- < y 5i (x-L 5] ) --^(x-L^ 
+—co, (x-L^ ) --co,, (x-L, \ +—co v (x-L, ) -—co^ (x-L, ) 

9 J 2 \ J 12 / O 3 \ °>12 / O 3 \ ^123 / 9 4 \ J 123 / 



(17d) 



Then, the deflection ysj(x) in each section of the multilayered microresonator is determined by 
substituting Equation (17) into Equation (9) and considering the boundary conditions expressed by 
Equation (10), as well as the integration rules of the Macaulay's functions. Thus, the deflections ys/x) 
can be determined as: 
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For 0<x<L„ 
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For L„ <x<L, 
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'1234 
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For multilayered microresonators with two or three different sections, only the first two or three 
deflections equations must be used. Next, ysj(x) is used into Equations (6) and (7) to find the maximum 
potential and kinetic energy of the microresonator. Finally, by substituting these energies into Equation (8), 
the first bending resonant frequency of the multilayered microresonator can be estimated. 
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3. Application of the Analytical Model 

We applied our analytical model to estimate the first bending resonant frequency and deflection of 
two multilayered microresonators reported in the literature. In addition, we used finite element models 
(FEMs) to study the bending vibration of these microresonators. Our analytical results agree with the 
FEMs and experimental results of the two microresonators. 

We considered the geometrical configuration of two multilayered microresonators for mass sensing 
applications developed by Lu et al. [46,47]. Figure 5 shows a SEM image of the first multilayered 
microresonator of Lu et al. [46], which we called microresonator type-A. In addition, we named their 
second multilayered microresonator microresonator type-B (see Figure 6). 

Figure 5. SEM micrograph of the multilayered microresonator type-A [46]. Reprinted with 
permission from Elsevier Science B.V. Copyright© 2009. 




Figure 6. SEM micrograph of the multilayered microresonator type-B [47]. Reprinted with 
permission from Japan Society of Applied Physics. Copyright© 2007. 




These microresonators have self-actuation and self-sensing capability. The details of their 
fabrication process can be found elsewhere [46,47]. Figure 7 shows the geometrical configuration of 
the layers and load types on the microresonator type-A. It uses small electrodes of piezoelectric lead 
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zirconate titanate (PZT) for microresonator excitation and contains a piezoresistive gauge for detection 
of its resonant frequency (see Figure 8). It can be divided in four different sections of layers. The first 
two sections contain nine layers composed by the following materials: silicon (Si), silicon dioxide 
(Si0 2 ), titanium (Ti), platinum (Pt), PZT, chromium (Cr), and gold (Au). The order of the nine layers 
on the first and second section is the following: Si/Si0 2 /Ti/Pt/PZT/Cr/Au/Cr/Si0 2 . Only the first 
section has two holes to decrease initial stress due to PZT deformation, which improve the signal to 
noise of the piezoresistive gauge contained in the first section. The third only contains a Si layer and 
the fourth section has three layers of Si/Cr/Au. 

Figure 7. Geometrical configuration of the layers and load types on the microresonator 
type-A. This figure is not drawn in scale. 




Figure 8. Detail view of the layers located on the first three sections of the microresonator 
type-A. This figure is not drawn in scale. 



Electrodes Position of piezoresistive gauge 




section 

Table 1 indicates the order of the layers in each one of the four sections for the microresonator 
type-A. Figure 9 depicts the main dimensions of the microresonator type-A. In addition, Table 2 
indicates the geometrical parameters of all layers used in the microresonator type-A. For the first 
section, bis\ include the overall width of the layers located to a same distance his\. 
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Table 1. Order of the layers in the four sections of the micro resonator type- A. 



Layer 


Section 1 


Section 2 


Section 3 


Section 4 


Si 


1st 


1st 


1st 


1st 


Si0 2 


2nd 


2nd 






Ti 


3rd 


3rd 






Pt 


4th 


4th 






PZT 


5th 


5th 






Cr 


6th 


6th 




2nd 


Au 


7th 


7th 




3rd 


Cr 


8th 


8th 






Si0 2 


9th 


9th 







Figure 9. Detail view of the layers located on the first three sections of the microresonator 
type-A. This figure is not drawn in scale. 



50 urn 100 jam 50 |am 500 (am 




Table 2. Geometrical parameters of the layers used in the microresonator type-A. 



Geometrical 
parameter 


Dimension 
(urn) 


Geometrical 
parameter 


Dimension 
(urn) 


Ls\ = Lsi 


50 


blSi = &3S4 


130 


Lsi 


100 


hsi = hsi = hsi = hsA 


5.0 


Lsa 


500 


hsi — hsi 


0.40 


bm 


140 


hsi ~ hsi 


0.05 


bis,\ 


120 


hsi = hsi 


0.25 


blS\ = ^4S1 = ^5S1 = ^6S1 


72 


hsi ~ hsi 


1.0 


blS\ = = = &3S2 


72 


hsi = hsi = hs\ = hsi 


0.03 


&4S2 = ^552 = ^652 = &7S2 


72 


hs\ = hsi = hsA 


0.15 


biSl = &9S2 


72 


hsi = hsi 


0.30 


b\Sl = &1S4 


150 







Figure 10 shows the geometrical configuration of the layers and load types on the microresonator 
type-B, which contains four different sections of layers. Its first section is integrated by nine layers that 
are fabricated with the same materials used in the first section of the microresonator type A. The 
second section is composed by two layers (Si/SiCh) and the third section is only formed by Si (see 
Figure 11). The fourth section contains three layers that are made of Si/Cr/Au. The order of the layers 
in each one of the four sections for the microresonator type-B is showed in Table 3. Table 4 indicates 
the geometrical parameters of all layers employed in the microresonator type-B. 
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Figure 10. Geometrical configuration of the layers and load types on the microresonator 
type-B. This figure is not drawn in scale. 




Figure 11. Detail view of the layers located on the first three sections of the 
microresonator type-B. This figure is not drawn in scale. 



9 layers. 




section 3rd 

section 

Table 3. Order of the layers in the four sections of the microresonator type-B. 



Layer 


Section 1 


Section 2 


Section 3 


Section 4 


Si 


1st 


1st 


1st 


1st 


Si0 2 


2nd 


2nd 






Ti 


3rd 








Pt 


4th 








PZT 


5th 








Cr 


6th 






2nd 


Au 


7th 






3rd 


Cr 


8th 








Si0 2 


9th 
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Table 4. Geometrical parameters of the layers used in the microresonator type-B. 



Geometrical 
parameter 


Dimension 
(um) 


Geometrical 
parameter 


Dimension 
(^im) 




80 


hm = h lS2 = h lS 3 = h m 


5.25 










j 


1 A 

14 


7 


0.05 


Lsa 


100 




0.25 


b\S\ = = b\si = ^252 


50 




1.0 


b?,S\ = &4S1 = = b(,S\ 


38 


hesi = ^851 = ^254 


0.03 


^751 = ^851 = bgsi 


38 


^751 


0.15 










blS3 = b\S4 


50 


^354 


0.24 



Both microresonator type-A and type-B have four sections with different numbers of layers in each 
one of them. In addition, their layers have a symmetrical configuration with respect to the xy-plane. 
Therefore, our analytical model for a general multilayered microresonator can be used to estimate the 
first bending resonant frequency of both microresonator type-A and type-B. Table 5 shows the 
mechanical properties of all layers used in the two microresonators, which were obtained from the 
literature [48-51]. Based on the mechanical properties and values of the geometrical parameters of the 
two microresonators, we first find the weight by unit length, reaction load, bending moment, and 
effective stiffness for both microresonators (see Table 6). Using these values, we substitute the 
deflections Equation (18) into Equations (6) and (7) in order to determine the maximum potential and 
kinetic energies. Finally, these energies were applied into Equation (8) to estimate the first bending 
resonant frequency (f r ) of the microresonators. Thus, our analytical results for the microresonator type-A 
and type-B are 13.18 kHz and 186.51 kHz, respectively. These results agree well with the experimental 
results of both microresonators obtained by Lu et al. [46,47], which are 12.81 kHz and 186.0 kHz, 
respectively. Our analytical results present a small relative difference of 2.9% and 0.3% with respect to 
experimental data. Also, we compared our analytical results with those obtained through two FEMs 
and a simple cantilever model (SCM). 

Table 5. Mechanical properties of the layers used in the microresonator type-A and 
microresonator type-B. 



Layer 


Young's modulus 


Density 


Poisson 






(kg-nT 3 ) 


ratio 


Si 


130 


2330 


0.28 


Si0 2 


73 


2200 


0.17 


Ti 


110 


4510 


0.32 










PZT 


63 


7500 


0.36 


Cr 




7190 


0.21 


Au 


75 


19,300 


0.42 
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Table 6. Values of the weight by unit length, reaction load, bending moment, and effective 
stiffness for the microresonator type- A and microresonator type-B. 



Parameter 


Microresonator 


Microresonator 




type-A 


type-B 


R 0 x 10 9 N 


15.90 


1.96 


M 0 x 10 12 N-m 




flfci x 10 6 N-m _1 


29.09 


12.79 


&fax K^N-m" 1 


30.32 


6.43 


&fox lO^-m" 1 


17.14 


6.00 


&fct x K^N-m" 1 


21.11 


8.09 


(EI z ) sl x 10 12 N-m 2 


344.98 


159.86 






(EI z )siX 10 12 N-m 2 


203.12 


78.38 









For a simple cantilever model where the origin of its coordinate system is located at its fixed end 
and the xz-plane is considered its neutral plane, the first bending resonant frequency (f rc ) can be 
calculated as [52]: 

f rc =0.162^15 (19) 

where h c and L c are the thickness and length of the cantilever, respectively; E c and p c are the Young's 
modulus and density of the cantilever material. 

The deflection curve (y c (x)) for a simple cantilever model can be determined as [53]: 

c c 

where co c and I c are the weight by unit length and second moment of area of the cantilever. 

Both microresonators type-A and type-B can be approximated to a simple cantilever model if only a 
layer in each one of their four sections is considered. For this, we used the layer with the largest 
thickness, which corresponds to the silicon layer for both microresonators. For the microresonator 
type-A, we neglected the holes and considered a silicon layer with dimensions 700 x 150 x 5 |im. This 
layer presents dimensions much larger than any other layer of the microresonator type-A. For the 
microresonator type-B, the dimensions (80 x 38 x 1 |im) of its piezoelectric layer (PZT) are close to 
those of its silicon layer (200 x 50 x 5.25 |im). Using Equation (19), the first bending resonant 
frequencies for both microresonators are 12.35 kHz and 158.82 kHz, respectively. For the 
microresonator type-A, the SCM's result has a small relative difference of -3.6% with respect to 
experimental result. For this case, the approximation of the microresonator to a SCM is suitable 
because the microresonator' s dimensions are close to dimensions of its silicon layer {i.e., the 
piezoelectric and isolation layers have dimensions much smaller than those of the silicon layer). 
Nevertheless, the SCM's result for the microresonator type-B has a high relative difference of -14.6% 
with respect to the measured frequency. This high relative difference is caused because both metallic 
and isolation layers have an overall thickness and length (2 |im and 80 |im) close to those of the silicon 
layer (5.25 |im and 200 |im). Therefore, the mass and stiffness of the metallic and isolation layers have 
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a higher influence on the mass and effective stiffness of the microresonator. For this case, our 
analytical model is more suitable than the SCM. 

In addition, we made two FEMs for both microresonators through ANSYS® software. These models 
used solid95 type elements, in which each element is defined by 20 nodes with three degrees of 
freedom per node: translations in the nodal x, y, and z directions. For FEM of the microresonator type-A, 
we only considered its thicker layers (silicon, SiCh, and PZT) in order to obtain a mesh that does not 
overcome the maximum nodes number (125,000 nodes) allowed by our ANSYS® software license. 

Figure 12. First bending mode of the microresonator type-A obtained using a FEM. 



AN 

NODAL SOLUTION 
STEP=1 




0 .222225 .444451 .666676 .888901 

.111113 .333338 .555563 .777789 1 



Figure 13. First bending mode of the microresonator type-B obtained using a FEM. 



AN 

NODAL SOLUTION 

STEP=1 
SUB =1 
FREQ=188416 




0 .222269 .444537 .666806 .889075 

.111134 .333403 .555672 .77794 1 



The thinner layers of the microresonator type-A have thicknesses less than 300 nm, which 
significantly increases the number of nodes in the FEM mesh. This is a problem to mesh layers on the 
order of nanometers. Similarly, the FEM of the microresonator type-B only included the thicker layers 
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(silicon, SiC>2, and PZT layers). Figures 12 and 13 show the first bending vibration mode of the FEMs, 
which were obtained to 13.49 kHz and 188.42 kHz. These results have a relative difference of 5.3% 
and 1.3% with respect to the experimental results. 

In addition, we compared the normalized deflection y(x)/y max of both microresonators obtained by 
the three models studied in this work. For this, we used the deflections Equations (18) and (20), 
dividing their results by the maximum deflection. Figures 14 and 15 illustrate the normalized 
deflections along the length of both microresonators, which were obtained through the SCM, FEMs, 
and our analytical model. The results of normalized deflections obtained by our analytical model agree 
very well with those of the FEMs. 

Figure 14. Normalized deflection j(x)/y max of the microresonator type-A obtained through 
our analytical model, a FEM, and a simple cantilever model. 




Length of microresonator type-A (|im) 

Figure 15. Normalized deflection y(x)/y max of the microresonator type-B obtained through 
our analytical model, a FEM, and a simple cantilever model. 




Length of microresonator type-B (jim) 
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Our analytical model can be used in the design phase of multilayered microresonators with variable 
cross section in order to estimate their lowest bending resonant frequencies and deflections. Using our 
model, a designer can determine the dimensions of the microresonator's layers that allow it to operate 
at a resonant frequency suitable to a particular application. A designer can also use the proposed model 
to know the influence of the microresonator's materials on its resonant frequency. 

4. Conclusions 

An analytical model to estimate the first bending resonant frequency and deflection curve of 
multilayered microresonator with variable cross-section was presented. This model is formulated 
through the Rayleigh and Macaulay methods, as well as the Euler-Bernoulli beam theory. We have 
applied our analytical model to two multilayered microresonators composed by layers of seven 
different materials reported in the literature. The results of proposed analytical model presented a 
relative difference less than 3.0% with respect to experimental data. Our analytical model can be useful 
in the mechanical design of multilayered microresonator for detection of mass and chemical species. 
Future work will include the effect of the residual stress on the multilayered microresonators with 
variable cross-section. In addition, we will study the resonant frequency shift of these microresonators 
due to films or particles deposited on their surface. 
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